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Abstract. In an attempt to determine the nature of the 
enigmatic cocoon stars in the Quintuplet Cluster, we have 
obtained mid-infrared imaging and spectrophotometry of 
the cluster, using the CAM and SWS instruments on ISO, 
using SpectroCam-10 on the Palomar 5m telescope, and 
NICMOS on HST. The spectra show smooth continua 
with various dust and ice absorption features. These fea- 
tures are all consistent with an interstellar origin, and 
there is no clear evidence for any circumstellar contribu- 
tion to these features. We find no spectral line or feature 
that could elucidate the nature of these sources. Detailed 
modeling of the silicate absorption features shows that 
they are best reproduced by the fj, Cep profile, which is 
typical of the interstellar medium, with t s ii ~ 2.9. The 
high spatial resolution mid-IR images show that three of 
the five cocoon stars have spatially extended and asym- 
metric envelopes, with diameters of ~ 20, 000 AUs. 

A reddening law similar to that of Lutz (1999) but with 
silicate features based on the fj, Cep profile and normalized 
to our value of t s h is used to deredden the observed spec- 
trophotometry. The dereddened energy distributions are 
characterised by temperatures of 750-925 K, somewhat 
cooler than determined from near IR data alone. Models 
of optically thin and geometrically thick dust shells, as 
used by Williams et al. (1987) for very dusty, late- type 
WC stars, reproduce the observed SEDs from 4 to 17/itm, 
and imply shell luminosities of log(L/L Q ) ~ 4.5-4.9 for 
the brightest four components. An analysis of the various 
suggestions proposed to explain the nature of the cocoon 
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stars reveals serious problems with all the hypotheses, and 
the nature of these sources remains an enigma. 
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1. Introduction 

The Quintuplet Cluster (AFGL 2004) is one of the three 
young clusters known in the vicinity of the Galactic Cen- 
tre. It was first identified in the near infrared survey of 
Glass, Catchpole, & Whitelock (1987) as a bright source 
coincident with AFGL 2004, and found to be a bright 
IRAS source with a cool energy distribution (Glass 1988). 
It is located about 30 pc from the Centre in projec- 
tion, and its cluster nature became apparent with the 
advent of IR cameras through the works of Glass, Mon- 
eti, & Moorwood (1990; GMM90), Nagata et al. (1990), 
Okuda et al. (1990), Moneti, Glass, & Moorwood (1992), 
and Moneti, Glass, & Moorwood (1994 [MGM94]), who 
identified over a dozen stars in the cluster, including 
five with extraordinarily large luminosities (~ 10 5 Lq), 
cool (600 — 1200 K) energy distribution, and featureless 
if-band spectra. These characteristics led those authors 
to tentatively identify them as massive, dust-enshrouded 
young stars, henceforth cocoon stars. The possibility that 
they be OH/IR stars was discarded since (i) they have 
no CO bandhead at 2.3 and 4.6 /xm, (ii) they show little 
or no variability (Glass et al. 1999 find peak-to-peak of 
~ 0.45 ± 0.05 mag and K for two of them, no variability 
for the remainder), (iii) they lack OH maser emission, (iv) 
they are cooler and more luminous than typical OH/IR 
stars, and (v) they are much younger than OH/IR stars, 
assuming they are coeval with the cluster (« 4 Myr, Figer 



2 



A. Moneti et al.: The Cocoon Stars in the Quintuplet Cluster 



McLean and Morris 1999 [FMM99] , Figer et al. 1999). 
More recently, FMM99 suggested that the cocoon stars 
might be extremely dusty late-type WC stars (DWCL, 
Williams et al. 1987, WHT), even though they do not de- 
tect the expected emission features in the J-band. And 
finally, Glass et al. (1999) suggested they could be self- 
obscured, very massive O stars close to the ZAMS, which, 
according to Bernasconi and Maeder (1996), go through 
a prolonged accretion phase of 2-2.5 Myr, comparable to 
the age of the Quintuplet Cluster, before becoming visible. 

In an attempt to elucidate the nature of the cocoon 
stars, we have obtained spectroscopy and spectrophotom- 
etry in the mid-IR, where these sources are brightest, in 
order to (i) search for spectral features intrinsic to the 
stars, (ii) obtain accurate measurements of their spectral 
energy distribution, (iii) to investigate their spatial exten- 
sion, and (iv) to study the shape of the silicate feature in 
particular to determine whether it could be partly intrin- 
sic to the sources. 

2. Observations, calibration, and results 

Two instruments onboard the Infrared Space Observatory 
(ISO, Kessler et al. 1996) were used: ISOCAM (Cesarsky 
et al. 1996) was used to obtain high spatial and low spec- 
tral resolution Circular Variable Filter (CVF) spectropho- 
tometry of the Quintuplet Cluster from 2.5 to 17/mi, and 
SWS (de Graauw et al. 1996) was used to obtain two high 
resolution spectra at two positions in the cluster. Table 1 
lists the ISO identification number (ION) of the datasets 
used. Of these, the last was proposed and planned by the 
authors, while the remainder were obtained by other in- 
vestigators and taken from the archive. Some preliminary 
results based on the ION 09901781 dataset have already 
presented Nagata et al. (1996). Since that time, the cal- 
ibrations have improved considerably and, in particular, 
the CVF spectral response function used here is based on 
stellar observations obtained in orbit, and differs consid- 
erably from the one used in Nagata et al. (1996), which 
was obtained on the ground prior to launch. 

We also present high spatial resolution ground-based 
8-14 /jm spectroscopy and narrow-band imaging photome- 
try of the cluster members, obtained with SpectroCam-10 
on the Palomar 5 m telescope, and 1-2.5 /mi photome- 
try derived from HST/NICMOS images. These data are 
particularly useful to understand the limits of low spatial 
resolution ISO data. 

2.1. SWS Grating Spectroscopy 

The SWS spectra were obtained in the AOTl-speed 3 
mode, which scans the grating continuously over its full 
range, and produces a number of spectral segments, cor- 
responding to the various aperture/grating/detector com- 
binations. These segments cover the full spectral range of 
2.3 — 47 /im with a nominal resolving power ranging be- 
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Fig. 1. Location of the SWS apertures on a 7.4 /tm ISO- 
CAM image of the Quintuplet Cluster. The intensity scale 
is logarithmic, north is up and east is to the left. The nam- 
ing convention is that of GMM90. Solid line is aperture 1 
(14" x 20"), and dashed lines is aperture 2 (14" x 27"). 
The source to the NE of Q4 was denoted Q5 in GMM90. 



Table 1. ISO Datasets 



ISO ION 


Instr. 


mode 


date 


09901781 


CAM 


C04 (CVF) 


24-Feb-96 


28701246 


SWS 


SOl-speed 3 


29-Aug-96 


29702147 


SWS 


SOl-speed 3 


09-Sep-96 


84500303 


CAM 


C04 (CVF) 


09-Mar-98 



tween 250 and 600. In this work we will be concerned 
primarily with data out to 28 /im, which was obtained 
with the two smallest apertures, which are nominally 
14" x 20", 14" x 27", respectively, though in practice 
slightly larger [A. Salama, private communication]. The 
aperture changes at 12 /im. The location of the nominal 
apertures on a 7.4 /im images CVF image is shown in Fig- 
ure 1, where the cocoon stars are identified following the 
naming convention of Moneti, Glass, & Moorwood (1992). 
The SWS apertures were centered on Q4 and Q3 (respec- 
tively GCS 3-1 and GCS 4 in Kobayashi et al. 1983). 

Data processed with version 7.0.1 of the ISO Pipeline 
were further cleaned with ISAPjj] by removing data from 
noisy detectors, and, for each spectral segment, the data 
from the different detectors were normalized to the aver- 
age of the complete set of data, effectively compensating 

1 The ISO Spectral Analysis Package (ISAP) is a joint de- 
velopment by the LWS and SWS Instrument Teams and Data 
Centers. Contributing institutes are CESR, IAS, IPAC, MPE, 
RAL and SRON. 
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Fig. 2. SWS spectra of Q2 (primarily) and Q3. The jump at 12 /im corresponds to the aperture change; see text for 
details. 



for slight residuals in dark current (at low signals) or in de- 
tector responsivity (at high signals). Following these steps, 
the agreement between adjacent segments was quite good; 
the results are shown in Figure 2. The jump at 12 /im is 
due to the change in effective aperture, and is expected 
when observing extended sources. At wavelengths beyond 
28 /im the effective SWS apertures are even larger and 
they include a significant part of the Pistol Hn region 
whose hot dust emission dominates the spectrum. These 
parts of the spectra will not be considered further. 

The calibration used in the SWS Pipeline assumes that 
the target is point-like and in the center of the aperture. 
In that case, the nominal overall uncertainties would be 
- 5% up to 4.1 ^m, ~ 10% for 4.1— 27 /im, and - 30% 
beyond 29 /im (see de Graauw et al. 1999 for details). This 
is essentially the case for Q3, which is sufficiently isolated 
from the other sources. In the case of the GCS-3 complex, 
which includes Ql, 2, 4, and 9, the aperture was centered 
on Q4, but the brightest source in the beam is Q2, located 
6" south of Q4, and significant light is scattered into the 
aperture from the other components, especially at long 
wavelengths, where the PSF is broadest. 

2.2. CAM CVF Imaging 

The ISOCAM instrument, used to obtain CVF (Contin- 
uously Variable Filter) scans, contains two separate cam- 
eras that cover the regions shortward and longward of 
5 /im (henceforth 'SW and 'LW, respectively). Both cam- 
eras utilize 32 x 32 pixel arrays, an InSb array with CID 



readout for the SW channel and a Si:Ga array with a di- 
rect read-out for the LW channel. The spectral resolving 
power ranges from 45 to 100. For each CVF segment, a 
single scan was obtained with wavelength increasing with 
time, and each scan was preceded by a long (40 — 100) 
series of frames obtained at the first CVF position for de- 
tector stabilisation. 

We present here 2.5-3.0 /im and 4-5 /im spectra from 
the ISO 09901781 dataset, and 5-17 /im spectra from the 
ISO 84500303 dataset. All these data were obtained with a 
pixel scale of 1.5". The ISO 09901781 dataset also contains 
5-17 /im data obtained with the 3 arcsec/pixel scale, but 
these data are, for the most part, heavily saturated, and 
are not discussed here. In both cases, the total integration 
time per spectral position was approximately 1 min. 

For completeness we also examined the ISO 49400231 
and the ISO 49400268 datasets, which contain 3x3 rasters 
of the cluster at two CVF positions and with two broad 
filters. The filter data suffer from significant saturation, 
and the CVF data are qualitatively of much lower quality 
than the complete CVF scans, probably as a results of 
transient effects during the raster scanning process. 

The ISOCAM data were reduced and calibrated with 
the IDL-based CAM Interactive Analysis (CIA) package^ 
following the standard steps of (i) dark subtraction, (ii) 



2 CIA is a joint development by the ESA Astrophysics Divi- 
sion and the ISOCAM consortium, led by the ISOCAM PI, 
C. Cesarsky, Direction des Sciences de la Matiere, C.E.A., 
France 
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Fig. 3. Composite of most of the CVF images of the Quintuplet Cluster. Wavelengths in /im are labeled. Individual 
panels are normalized in XF(X) units. North is up and east is to the left, and each panel is 48" on the side. The jump 
in image position between 5.12 and 5.26 /xm is due to the switch from the SW to the LW section of ISOCAM. 
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removal of cosmic ray events, (iii) correction of the be- 
haviour of the detector in response to signal jumps (tran- 
sient correction, LW only, using the "inversion" method), 
(iv) averaging all valid images obtained at a given CVF 
position, (v) flatficlding the results using library flats, and 
(vi) conversion the average signal to physical units using 
standard conversion factors, originally derived from sev- 
eral standard stars. Throughout the analysis, the dead 
column (#24) in the LW array was replaced by the av- 
erage of the adjacent columns. No ghost image could be 
identified in the CVF data, nevertheless some light could 
be lost to low-level, diffuse ghost images. From a study of 
the background, sampled in the NE corner of the array, 
see §3.2, we find that the background level is 2% of the 
total signal in the array at A < 6 /xm, then rises quickly to 
5% at 8 /xm, and remains at that level thereafter. 

A composite of the resulting data is shown in Figure 
3, where the individual panels are normalized in XF(X) 
units. Note in particular (i) the nearly complete silicate 
absorption, (ii) that the brightest source at short wave- 
length, Q2, is no longer so at long wavelength, where Q9 
becomes brightest, and (iii) the increase in PSF size with 
wavelength. 

The nebula south of the cluster, which becomes promi- 
nent at A > 10 /im, is the Pistol nebula, and the Pistol 
star is prominent at the bottom of the frames at the short 
wavelengths. These are the subject of a separate paper 
(Moneti ct al. 2000, in preparation). 

The main difference between these results and those of 
Nagata et al. (1996) is that the spectral response function 
for the SW-CVF section (A < 4 /an) has changed consid- 
erably (Biviano et al. 1997), so that the new flux densities 
in the 2.5-3 /<m range are now lower by almost a factor 
of 2, and the 3 /im ice feature is shallower than previously 
reported. As will be shown below, the new results are in 
better agreement with previously published results. 

2.3. Palomar Imaging and Spectroscopy 

High spatial resolution (1") images and low-resolution 
(R ~ 50) spectra were obtained in the 8-13 /xm at- 
mospheric window with SpectroCam-10 (Hay ward et 
al. 1993) on the Palomar 5m Hale telescope. These ob- 
servations were taken on July 10-12, 1993 (spectroscopy 
and 11.7 /xm imaging) and June 14, 1995 (8.7 /xm imaging). 
SpectroCam-10 uses a 128 x 128 Rockwell Si:As detector. 
Its camera mode has a fully illuminated field of view of 
~ 15" diameter with 0.256" pixels. All observations em- 
ployed the standard mid-infrared techniques of chopping 
the secondary mirror (typically at 5-10 Hz) and nodding 
the telescope (typically 60" E-W, every few minutes) in 
order to adequately subtract the thermal background. Im- 
ages were obtained with two standard OCLI filters cen- 
tered at 8.74 and 11.66 /xm with half-power bandwidths 
of AA = 0.78 and 1.11 /xm, respectively, and mosaics were 
constructed from these images. The flux calibrators were 
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Fig. 4. Mosaic of the Quintuplet Cluster and the Pis- 
tol Nebula at 11.7 /xm. Pixels are 0.256" wide, the grey- 
scale is logarithmic to show the large dynamic range, 
and the contour levels are 1.3, 3.1, 8.0, 120 49, 123, and 
311 x 10- 14 Wm- 2 /mi- 1 arcsec" 2 (0.05, 0.12, 0.31, 0.76, 
1.9, 4.8, and 12.1 Jyarcsec~ 2 ). 

a Boo and a Her, and their fluxes were obtained from 
Hanner et al. (1984). An airmass correction was derived 
from observations of standard stars at different airmasses; 
the cocoon stars were observed at a typical airmass of 2.3. 

Figure 4 represents a mosaic of 47 seeing-selected, flat- 
fielded 11.7 /xm images, each with a total on-source inte- 
gration time of 40 seconds. During the mosaic process the 
pixels were subsamplcd by factor of three in order to opti- 
mise image registration; after mosaicking, the images were 
rebinned back to the original pixel sizes, and the result 
was smoothed slightly with a Gaussian of a = 0.8 pixel. 
A mask was used to omit pixels not fully illuminated. The 
image registration accuracy in the region containing the 
cocoon stars is estimated to be better than 0.3 pixies, or 
< 0.1"; in regions without bright stars, the telescope off- 
sets were used, and the registration accuracy is only ~ 1". 
The 8.7 /xm image was constructed from 56 images with 5 
seconds on-source integration time each. 

The spectroscopy was obtained with a 2" x 15" slit to 
yield a resolution of A/AA ~ 50. The slit was oriented 
N-S, the dispersion was 0.046 /xm/pixel, and four pixels 
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were combined for each resolution element. Four grating 
settings were used to cover to 8-13 jum range and each set- 
ting was exposed for a total integration time on-source of 
20 seconds for each of the cocoon stars. Atmospheric lines 
were used for wavelength calibration and f3 Peg was ob- 
served for flux calibration and to correct for atmospheric 
lines. 

2.4. HST/NICMOS Imaging 

The Quintuplet cluster was observed with NICMOS on 
HST in program 7364 in Sept. 1997 with Camera 2, giv- 
ing a spatial resolution of 0.075"/pix. Three broadband 
filters were used: F110W, F160W, and F205W. Figer et 
al. (1999c) give details of the observations. The raw data 
were retrieved from the STScI archive and were re-reduced 
with an improved linearity correction algorithm using the 
'calnica' task in STSDAS. This had the effect of recov- 
ering photometric information for the central pixels of 
the Quintuplet stars, which were flagged as saturated in 
the pipeline data reduction. Mosaics were constructed by 
aligning 16 dithers in each filter to produce the final im- 
ages. Residual dc levels were measured and subtracted 
from the images. 

2.5. Photometric and Spectrophotometric Results 

Photometry of the cocoon stars was extracted from each 
CVF image by integrating within a small aperture and 
multiplying the result by an aperture correction fac- 
tor that was determined from CVF observations of the 
standard star 5 Draconis. The sky background at each 
wavelength was determined from regions of the image 
free of stars, and the same sky level was used for all 
sources. Even with a small aperture the photometry of 
the faintest sources can suffer significant contamination 
from the wings of the brighter sources nearby. This ef- 
fect is most serious at long wavelength where the PSF is 
largest. This problem is most acute for Q4, which is the 
weakest source and which is surrounded by three bright 
ones. 

The SpectroCam-10 photometry was extracted with 
a 5-7" (diameter) aperture, depending on the extended- 
ness of the source. The NICMOS photometry was derived 
from the energy enclosed in a 3-pixel radius (0.45" diam- 
eter) aperture, and then corrected to an infinite aperture 
using measured isolated stars to derive the aperture cor- 
rection. These correction factors were 1.31, 1.50, and 1.71 
for the F110W, F160W, and F205W filters, respectively. 
The results are listed in Table 2. We assign these values 
conservative uncertainties of 0.2 mag, which are due pri- 
marily to the photometric extraction method and to errors 
in the estimate of the residual dc levels. The actual uncer- 
tainty for the brighter sources is likely to be closer to 0.1 
mag, especially at the longer wavelengths. 



Conversion to flux densities are based on the system 
in which Vega (without its mid-IR excess) is assigned 0- 
mag m all bands (Cohen et al. 1992)0 While conversion 
of magnitudes to flux densities is straightforward, the de- 
termination of the proper wavelength at which the flux 
density should be plotted is more complicated, since it 
depends on the shape of the (observed) energy distribu- 
tion of the source. For this purpose, we have modeled the 
intrinsic shape of the 1-4 /jm energy distributions as warm 
blackbodies, and reddened them with the reddening law 
discussed below to produce a model observed energy dis- 
tribution. This was then used to compute the isophotal 
wavelength, Aj SO , as defined by the relation 

^(Aiso) J T(X)d\ = J F(X)T(X)d\, 

where F(X) is the observed spectrum of the source, and 
T(A) is the transmission of the instrument (including fil- 
ter, detector efficiency, optics, and possibly the earth's 
atmosphere). For Ai SO to be unique, F(X) must be a 
monotonic function of A, which is indeed the case for 
our models. The isophotal wavelengths for the NICMOS 
and SpectroCam-10 filters are listed in Table 2 together 
with the values of reddening correction {Ax/ Ay) and the 
adopted flux density for a zero magnitude FZM^/Jy. 

Figure 5 shows the observational results, which are 
compared with the ground-based photometry of FMM96 
(filled circles) and of Okuda et al. (1990, filled triangles, 
see also corrections in Nagata et al. 1996). Approximate 
Ai SO values were computed for the FMM99 data using stan- 
dard J, H, and K filter profiles, but not for the mid-IR 
data, which was obtained with narrow-band filters. 

The agreement between the CVF and the Palomar 
spectra is excellent, as is the agreement between the CVF 
spectra and the ground-based photometry. The discrep- 
ancy between the CVF and the Palomar spectra of Q4 at 
A > 10 /im is due to the contamination of the CVF spec- 
trophotometry of Q4 by its bright neighbours (see above). 
We consider the CVF spectrum unreliable there, and will 
instead rely on the Palomar results. 

At NIR wavelengths, the NICMOS fluxes fall below the 
ground-based results since the former can better isolate 
the target sources from other nearby sources (see Figure 
2 of Figer et al. 1999c) 

2.6. Detailed comparison of SWS and CVF spectra 

A comparison of the CVF and the SWS spectra is pre- 
sented in Figure 6. The SWS spectra were transformed 
to surface brightness using the measured aperture sizes, 
which are 1.12x and 1.29x, and 1.06x larger than the 
nominal apertures for Bands 1,2, and 3, respectively (A. 

3 The Vega model used for the absolute calibration can be 
found in www.iso.vilspa.esa.es/users/expl_lib/ISO/wwwcal/- 
isoprep/cohen/vega9605.dat. 
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Fig. 5. Observed spectra and photometry of the cocoon stars. The J-band data and some H-band data were left off 
scale in order to avoid overly compressing the ordinate, which would make the spectral features too small to see (these 
data points are plotted in Figure 11 below). The CVF spectra are the blue histogram, Palomar spectra and photometry 
are orange histogram and squares, NICMOS photometry are green asterisks. The green circles are near-IR imaging 
photometry from FMM99, and the green triangles are mid-IR aperture photometry from Okuda et al. (1990). The 
panels show the sources roughly in their positions on the sky. 

Table 2. Photometry (in magnitudes) of the cocoon stars 



Source 


CCS 


[F110W] 


[F160W] 


[F205W] 


[8.7] 


[11.7] 


Ql 


3-IV 


16.63 


11.48 


8.37 


1.14 


0.13 


Q2 


3-II 


14.46 


9.71 


6.80 


0.54 


-0.38 


Q3 


4 


15.31 


10.62 


7.60 


1.18 


0.49 


Q4 


3-1 


15.78 


11.11 


8.35 


2.61 


1.94 


Q9 


3-III 


18.88 


13.54 


9.86 


0.80 


-0.46 


A iao / 




1.14 


1.63 


2.10 


8.69 


11.55 


Ax/Av 




0.285 


0.177 


0.119 


0.073 


0.081 


FZMa/Jv 




1571 


1075 


695 


50.4 


28.9 



Salama, private communication). These factors, derived 
from beam profiles, are 0th order approximations of a pa- 
rameter that is a function of wavelength. The change of 
correction factor introduces a ~ 10% jump at the Band 



1 /Band 2 interface (4 /im), where the nominal aperture re- 
mains unchanged. The CAM spectra were extracted using 
a synthetic rectangular apertures approximating the mea- 
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Fig. 6. Comparison of SWS spectra (bluw histogram) with spectra extracted from CVF images (red and orange crosses, 
indicating different CVF sections). 



sured apertures of SWS, and centered as shown in Figure 
1. 

Over the 4-12 /im range, SWS Band 2, the agreement 
between the SWS and the CVF spectra is excellent for 
Q3, which is fairly isolated. The SWS spectrum of Q2 
falls short of the CVF result; this difference is due to the 
fact that Q2 is near the edge of the SWS beam, where the 
responsivity is much lower than in the beam centre. 

For Band 3, the SWS shows an increase in surface 
brightness for both spectra. No such increase is expected 
for Q3 on the basis of the CVF spectra, while a minor one 
is expected for Q2. The difference is likely to come from 
the detailed shape of the SWS beam profile, which is not 
accounted for in this simulation. 

In Band 1 the CVF spectra show a rapid decrease with 
wavelength, which is unexpected and contrary to the SWS 
behaviour. We note that this region corresponds the the 
beginning of the CVF spectrum, where detector transient 
effects are most serious (no transient correction was ap- 
plied to the SW data since the behaviour of these detectors 
was not sufficiently well characterised), and to the edge of 
the SW-CVF, where the determination of the response 
function is most uncertain. We consider this portion of 
the CVF spectra unreliable. 

We conclude that, with the exception of the first part 
of the CVF spectrum, the agreement between the two 
datasets is excellent, and this suggests that the overall 



10% over most of the 



uncertainties could be less that 
spectral range. 

3. Results 

3.1. Emission lines 



One of the objectives of the SWS spectra of the cocoon 
stars was to search for emission lines that could clarify 
the nature of these stars. A careful search of the spectra 
has revealed only ionized hydrogen and forbidden lines at- 
tributable to the Pistol H II region. These are stronger in 
the spectrum of Q3, which is closer to the Hn region so 
that more line (and continuum) emission is included in the 
SWS aperture. In contrast, only the strongest lines are de- 
tected in the Q2 spectrum, and these are detected mostly 
at the longest wavelengths, where the SWS aperture is 
largest. Furthermore, these lines are not detected in the 
CVF and Palomar spectra, albeit at lower resolution, con- 
firming that they are not intrinsic to the cocoon stars. In 
contrast, these lines are clearly detected in the CVF and 
Palomar spectra of the Pistol H n region. Beyond 28 /im 
the spectra (not shown) were obtained with even larger 
apertures, and the contamination by the Pistol Nebula, 
both in terms of forbidden line emission and of thermal 
emission from the warm dust (Moneti et al. 2000), begins 
to dominate over the cocoon stars. Thus, no emission line 
attributable to the cocoon stars themselves was detected. 
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3.2. Dust and Gas Features 

The spectra contain many dust absorption bands; these 
have been the subject of a number of recent studies 
(Schutte et al. 1998, Gerakines et al. 1999, Chiar et 
al. 2000), while Moneti and Cernicharo (2000) have stud- 
ied the gaseous components (all these results come from 
the same SWS spectra used here) . These results show that 
the lines of sight toward Q2 and Q3 contain (i) cold molec- 
ular cloud material, traced by absorption by cold (~ 10 K) 
CO and H 2 gas and by various ice species (H 2 0, CO2, 
CO, and CH4), and more tenuous ISM material, traced 
by aliphatic and aromatic hydrocarbon features. Special 
attention has been paid to the narrow absorption feature 
at 6.2 /im, attributed to the C-C stretch in aromatic hy- 
drocarbons. Schutte ct al. (1998) suggested that it could 
be the UIR feature seen in absorption, and justify the 
lack of absorption at 3.28 and 11.3 /J,m on excitation and 
temperature grounds, while they find that the observed 
spectrum is consistent with the presence of some absorp- 
tion at 7.7 /im. The best detection of this absorption fea- 
ture is toward WR118, a DWCL star whose absorption 
spectrum is believed to due to diffuse ISM material only, 
and marginal detections are reported towards three other 
DWCL stars (Schutte et al. 1998). Chiar et al. (2000), 
noting the 6.2 ^m absorption in the cocoon stars, which 
might also be DWCL stars, suggested that this absorp- 
tion feature could be associated with the environment of 
the DWCL stars rather than to the diffuse ISM. However, 
Chiar et al. (2000) also note that WC stars are not likely 
to be significant contributors to the carrier (s) of the C- 
H stretch (3.28 yum band), so that if the 3.28 and 6.2 /im 
bands (and the other UIR bands) are indeed connected, 
it is unlikely that they arise primarily in the surroundings 
of late-type WC stars. 

In contrast to the cocoon stars, the CVF data clearly 
show the UIR bands in emission in the spectrum of the 
background, as shown in Figure 7. The intensity of this 
emission is only ~ 2 and 3% of the continuum of Q2 and 
Q3, respectively, which is weak compared to the measured 
absorption. The background emission also shows the for- 
bidden lines of [Aril] and [Neil], and an important level 
of continuum. A part of this emission could be due to light 
scattered within the instrument. 

We have attempted to locate the source of the 6.2 /xm 
emission from the CVF data by computing the difference 
of the line planes of the datacube minus the adjacent con- 
tinuum planes. The result was flat, indicating that the 
emission is uniform over the whole region (unlike the [Ar 11] 
and [Ne 11] emission, which occur in the Pistol H 11 region). 
We deduce that the UIR emission is a characteristic of the 
line of sight rather than of the cocoon sources. We cannot 
determine where along the line of sight the UIR emission 
originate on the basis of our data. 
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Fig. 7. The UIR features as seen in the spectrum of the 
background, in this case a 5 x 5 pixel region in the north- 
east corner of the array. The UIR features at 6.2, 7.7, 8.6, 
and 11.3/um are identified by the arrows. 

3.2.1. The Silicate Features 

The silicate absorption bands at 9.7 and 18 fim are the 
strongest features in the cocoon star spectra. We have 
used several published cmissivity profiles to model the 
observed absorption. Most published profiles concern the 
9.7 fim feature which is observable from the ground: Roche 
and Aitken (1984, RA84), Hanncr, Brooke and Tokunaga 
(1995, HBT95), Bowey, Adamson, and Whittet (1998, 
BAW98), while two determinations include both features: 
Pegourie and Papoular (1995) and Simpson (1991). The 
studies of RA84 and BAW98 are based on absorption of 
ISM dust, HBT95 study the absorption by molecular cloud 
dust, and the remainder are based on the dust emission 
around late type supergiant stars, which are generally be- 
lieved to reproduce the ISM profile. These profiles, nor- 
malized to unity at their peak value, are shown in Figure 
8. A detailed comparison of these emissivities is beyond 
the scope of this work, but it is worth noting that there 
are important differences, even among those that should 
be very similar (e.g. the RA84 fi Cep profile and the Simp- 
son 1991 group 1 profile, which is based on several stars 
with narrow 9.7 /im features, including [i Cep itself). Some 
of the differences should probably be ascribed to the differ- 
ent methods of determining and modeling the continuum, 
which will produce important differences in the wings of 
the profiles. 

We removed the silicate features from the spectra by 
using the above emissivity profiles, e(A), as follows: 



F(X) = f(X)e T ^ (x '> 
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Fig. 8. Comparison of different published silicate emissivities. The adopted emissivity is traced by the black diamonds. 



where F(X) is the corrected spectrum, /(A) is the observed 
spectrum, and T a a is the peak optical depth of the fea- 
ture. The best fitting profile and the peak optical depth 
were determined subjectively by requiring that the cleaned 
spectrum be continuous and free of residual features. The 
alternative of using \ 2 methods was considered, but the 
errors at the beginning and the end of the profiles produce 
large deviations in the corrected spectra, which hide the 
deviations caused by more important errors in the profile 
where the emissivity is large. This is particularly evident 
over the 9.7 /im feature, which is strong and narrow. 

The best emissivity profile for the 9.7 /im feature is 
the /x Cep profile from RA84, and proper correction for 
all the sources is achieved with t s u = 2.9 ± 0.1. The weak 
and sharp feature remaining at 9.4 /im is not real, but it 
is due to a slight error in the determination of the emis- 
sivity profile in a region of rapid variation and where the 
atmospheric absorption is greatest. A change of 0.1 in t s h 
produces a clearly noticeable effect in the corrected spec- 
tra. 

The two 18 /im emissivity profiles are similar in shape 
but differ in strength relative to the 9.7 /im feature: 
Pegourie and Papoular (1985) determine r(18)/r(9.7) = 
0.35, while Simpson (1991) constrained her continuum 
model to produce t(18)/t(9.7) = 0.40. We have used the 
Pegourie and Papoular (1985) profile multiplied by 1.15 to 
produce r(18)/r(9.7) = 0.40, and joined it to the RA84 
yU Cep profile at 13 /im. We note that the Pegourie and 
Papoular (1985) profile of the 9.7 /im feature is much nar- 



rower than the cocoon star feature. This combined profile 
gives a deeper minimum at 12 /jm than if we had used the 
Simpson (1991) profile alone. The corrected (and dered- 
dened, see below) spectra, are shown in Figure 8. 

Our value of t s h is significantly higher than the value 
deduced by Schutte et al. (1998), who measured the depth 
of the feature relative to the level of the spectrum adjacent 
the feature. While their method is fine for their purpose, 
which is to intercompare the variation of this and other 
features in several sources, their result does not account 
for the absorption at 12-14 /im, as seen in Figure 8. 

We are confident that our method is reliable for the 
9.7 (im feature, but attempts to use to SWS spectra to 
place independent constraints on the strength and shape 
of the 18 /xm feature using the same criteria adopted for 
the 9.7 /im feature proved unsuccessful. The reason is 
probably that beyond 18-20 /im these spectra are a blend 
of various cocoon stars and of the Pistol H II region, and 
are thus intrinsically too complicated for a subjective de- 
termination of a "smooth" spectrum. Furthermore, the 
task was intrinsically more difficult due to the weaker and 
broader nature of the 18 /im feature. 

3.3. The Extinction Law 

The extinction law, A\/Ay, toward the GC is known to 
give rise to more extinction than the standard law both in 
the silicate bands and in the 5-8 /im spectral range: Roche 
& Aitken (1985, RA85) have shown that A v /T si i ~ 8.5, or 
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Fig. 9. Adopted extinction law compared to Draine & Lee (1984), Rieke, Rieke, & Paul (1989), and Lutz (1999). 



Ak/t s h ~ 0.70 with their extinction law (van de Hulst 
no. 15), toward the GC, which is about half the value in 
the solar neighborhood; and Rieke, Rieke, & Paul (1989) 
have shown on the basis of the Br-7 to Br-a ratios, that 
the extinction at 4 /im towards the GC is higher than ex- 
pected with standard laws. This results was confirmed and 
extended by Lutz (1999) using hydrogen emission lines 
throughout the 2.5 to 12 /im range. 

We have combined the above results to construct a 
reddening law appropriate for the GC as follows: 

1. shortward of the 2.5 /jm we follow the Rieke & Lebof- 
sky (1985) law which was determined for GC sources. 
At 0.9 /im < A < 2.5 /im this law can be approximated 
by a power law A\ ~ X~ a , with a = 1.63; 

2. between 2.5 and 3.8 /im some structure is seen in the 
Lutz (1999) results which is due to the ice feature at 
2.9 /im and to the amorphous carbon feature at 3.4 /im. 
We modeled these features as Gaussian; 

3. at A ~ 3.8-8 /im, a power law with a = 0.35, was used 
to approximate the results of Lutz (1999). Note that 
the point at 6 /xm, where the deviation from the power 
law is largest, suffers extra extinction due to the water- 
ice feature. At the position of Br-a this law gives the 
same results as the Rieke, Rieke, & Paul (1989); 

4. between 8 and 24 /im, the silicate absorption features 
are modeled using the emissivity adopted in the pre- 
vious section, with r s u = 2.9; and 



5. beyond 24 /im a power law with a = 2 is adopted 
(e.g. Draine & Lee 1984). In practice this last portion 
of the reddening law does not affect our results. 

The adopted extinction law is shown in Figure 9, where 
other commonly used reddening laws are shown for com- 
parison, and it is tabulated in Table 3. With this law, 
A V /A K = 8.9, A v /T s n = 10, and A K /r si i = 1.14. 

3.4- Extension of the Cocoon Stars 

Inspection of the Palomar images (Figure 4) shows signif- 
icant extended emission around some of the cocoon stars. 
Figure 10 shows the radial profiles which were measured 
on a special mosaic composed from selected exposures 
of superior seeing obtained during a single night at each 
wavelength. The radial profile of the isolated, but rather 
faint star, X Sgr, that was observed at nearly the same air- 
mass, is also shown for comparison (the profile of a Her 
observed at unit airmass, not shown, is even narrower, 
with a FWHM = 0.69" vs. 0.59" for the diffraction size of 
1.22A/D). 

The profiles of Q3 and of Q4 (out to 1.5", beyond which 
Q4 is likely to be contaminated by the bright sources 
nearby) are similar to that of reference star. Ql and Q9 
show wings that are almost 10 x higher, while Q2 shows a 
weaker, but clear excess. The profiles shown are azimuthal 
averages, but the extensions are clearly asymmetric and do 
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Table 3. The extinction law 



A/ /im 


A x /A v 


A/ um 


Ax/Av 


0.55 


1.0000 


12.00 


0.0721 


0.65 


0.7689 


12.50 


0.0634 


0.80 


0.5584 


13.00 


0.0569 


1.00 


0.4097 


13.50 


0.0545 


1.25 


0.2804 


14.00 


0.0553 


1.65 


0.1736 


14.50 


0.0561 


2.20 


0.1108 


15.00 


0.0573 


2.50 


0.0964 


15.50 


0.0589 


3.00 


0.0923 


16.00 


0.0612 


3.50 


0.0715 


16.50 


0.0643 


4.00 


0.0530 


17.00 


0.0677 


4.50 


0.0510 


17.50 


0.0707 


5.00 


0.0492 


18.00 


0.0726 


5.50 


0.0477 


18.50 


0.0730 


6.00 


0.0464 


19.00 


0.0721 


6.50 


0.0452 


19.50 


0.0704 


7.00 


0.0441 


20.00 


0.0683 


7.50 


0.0431 


20.50 


0.0663 


8.00 


0.0463 


21.00 


0.0644 


8.50 


0.0628 


21.50 


0.0626 


9.00 


0.0989 


22.00 


0.0610 


9.50 


0.1454 


22.50 


0.0586 


10.00 


0.1431 


23.00 


0.0561 


10.50 


0.1250 


23.50 


0.0537 


11.00 


0.1023 


24.00 


0.0515 


11.50 


0.0835 


> 24.0 


oc A" 2 



not show a common preferred direction, as seen in Figure 
4. 



Could the observed profile wings be due to other 
nearby red stars? The high spatial resolution 2.05 ^m 
(F205W) HST/NICMOS images of Figer et al. (1999) 
show three reddish stars close to Ql, a large blob of dif- 
fuse emission NW of Q2, weak but not red stars near Q3 
and Q4, and a fairly isolated Q9. The blob near Q2, the 
brightest of the group, is likely to be an artifact. We have 
performed photometry of all the "companions", and find 
that the red ones are far less red than the cocoon stars, 
and simple extrapolations of their energy distributions to 
8 and 12 ^im would produce negligible flux at those wave- 
lengths. Also, we note that Q9, which is the most extended 
cocoon star, is also the most isolated in the near IR. We 
therefore conclude that Ql, Q9, and to a lesser extend Q2 
are intrinsically extended. At the distance of the Galac- 
tic Centre (8 kpc), an angular radius of 2" corresponds 
to a physical radius of 0.08 pc or 16,000 Alls (radius to 
zero power). The true radius should be somewhat smaller, 
given the seeing of ~ 1" for the observations, but will 
remain nevertheless very important. 

3.5. Energy Distributions and Luminosities 

Using the reddening law constructed above we dereddened 
our spectra assuming Ak = 3.3 mag (Ay = 29 mag), as 
determined by Figer et al. (1998) from the average color 
excesses of the early-type stars in the cluster. The un- 
certainty in Ay is 5 mag. The results are shown in Fig- 
ure 11 together with the dereddened ground-based data 
of FMM99 and of Okuda et al. (1990) 

The resulting SEDs were modeled in three ways (see 
also Figure 11): 
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Fig. 11. Dereddened spectral energy distributions of cocoon stars; symbols are as in Figure 5. The models shown are 
simple blackbodies (dot-dash lines) and DWCL-like shells (solid lines). 



1. a constant temperature, optically thick shell, modeled 
by a Planck function (dot-dash lines in Figure 11); 

2. a constant temperature, optically thin shell, modeled 
by a Planck function multiplied by a A -1 emissiv- 
ity law: this produces an energy distribution narrower 
than the optically thick models (not shown); 

3. an optically thin but geometrically thick dust shell 
with r~ 2 density law (constant velocity wind) as used 
by WHT to model known DWCL stars. The shape of 
the model SED is determined by the dust tempera- 
ture, To, at the inner radius of the shell, i?o, and the 
shell thickness, H, which is measured in units of Rq. 
A dust emissivity proportional to A -1 was used. This 
model produces an SED that is wider than a Planck 
function modified with the same emissivity law, and 
whose width increases as the physical thickness in- 
creases (solid line in Figure 11). 

More detailed models of the DWCL dust shells have been 
computed by Zubko et al. (1998). While those models 
provide useful new results on the properties of the dust, 
the SEDs they produce are not significantly different from 



Table 4. Model parameters and results 



Source 




To/K 


H 


log(L/L ) 


Ql 


500 


550 


100 


4.5 


Q2 


625 


675 


50 


4.9 


Q3 


650 


700 


50 


4.6 


Q-i 


725 


775 


50 


4.1 


Q9 


400 


425 


100 


4.7 



those of WHT, so we will stay with the WHT models for 
the remainder of this work. 

The models described above were normalized to the 
data at A between 6.5 and 7 /mi, where the spectrum is 
reasonably free of dust or ice absorption features. The 
model parameters are listed in Table 4. In all cases the 
shell thickness is very large, and increasing H further does 
not have a significant effect on the SED, thus these models 
only place a lower limit to the shell thickness. A change of 
20 K in To produces a significant change in the SED, and 
the value of 20 K can be taken as the uncertainty in Tq. 
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As can be seen in the Figure 11, the DWCL-like mod- 
els can reproduce quite well the SEDs over the 4-17 /xm 
range. The models reproduce best the SEDs for Q2, which 
is the brightest in the mid-IR, and Q3, which is the most 
isolated source. Ql is not very bright at long A, and its 
spectrophotometry is likely to suffer some contamination 
from the nearby Q2, especially from the first Airy ring of 
the PSF of the latter. Q4 also suffers heavy contamina- 
tion by the nearby sources, as described in §2.5. Q9 is the 
the reddest source, its SED is somewhat peculiar, possibly 
due to some contamination at short A; for this source the 
model is mediocre at best. 

Below ~ 4 fxm the model SEDs fall far below the obser- 
vations, suggesting that the actual hot dust contribution 
is larger than in the model. Luminosities were computed 
integrating underneath the model SEDs; they are 10 4 ' 1 to 
10 4 ' 9 L Q for the DWCL-like models, while the blackbody 
models give luminosities about 0.1 dex lower (see Table 
4). These values do not differ significantly from previous 
determinations . 

4. Discussion 

4-1- The Origin of the Absorption Features 

Combining the adopted value of Ay with the value of 
T s ii = 2.9 ± 0.1 determined in §3.2, gives Ay /t b u ~ 10 ± 2. 
This, within the uncertainties, is equal to the nominal 
value of 8.5 determined by RA85 for the sources near Sgr 
A*. This implies that all the silicate absorption, and by 
extension all the extinction, takes place along the line of 
sight. This is not surprising: the infrared polarisation to- 
ward the cocoon stars and toward Sgr A* is about the 
same both in amount and in position angle (Kobayashi 
et al. 1983), and the depths of the other ice and gas fea- 
tures are all similar to, though usually slightly lower than, 
their corresponding values toward Sgr A* (Gerakines et 
al. 1999, Chiar et al. 2000, Moneti & Cernicharo 2000). 
The general conclusion is that there is slightly less absorb- 
ing material toward the Quintuplet Cluster than toward 
Sgr A* , and also that the fraction of cold molecular cloud 
material is lower toward the cluster. All this speaks against 
any absorption being produced by the cocoon stars them- 
selves, and supports the conclusion that all the absorption 
features observed are of interstellar origin. 

The lack of intrinsic absorption speaks against the 
interpretation of the cocoon stars as deeply embedded 
YSOs, where the deep silicate absorption is ascribed to 
the molecular cloud in which the YSO is deeply embed- 
ded. Furthermore, the 9.7 /xm feature toward deeply cm- 
bedded YSOs is usually broader than that of the ISM, and 
is better reproduced by the Trapezium profile than by the 
H Cep profile. Also, the spectra of embedded YSOs usually 
show warm CO and H 2 gas in absorption (see, e.g. van 
Dishoeck & Helmich 1996), presumably arising in a region 
where this gas is warmed either by the young stars itself 



or by shocks produced when the wind and/or the outflow 
from a young stars strikes the surrounding molecular ma- 
terial. No such warm gas is found toward the cocoon stars. 
Indeed, molecular line studies have shown the Quintuplet 
Cluster is located in a region free of molecular gas (Ser- 
abyn and Giisten 1991). 

4.2. On the DWCL Hypothesis 

The luminosities derived in §3.5 are those of the shells. 
The total stellar luminosity, assuming that the shell ab- 
sorbs and emits isotropically, will be the sum of the ob- 
served (and dereddened) shell luminosity and the observed 
(and dereddened) stellar luminosity, the latter being that 
fraction of the stellar luminosity not absorbed by the shell, 
and transmitted toward the observer. In the following we 
will alway refer to these observational quantities. 

For most of the DWCL stars studied by WHT, the ra- 
tio of shell to stellar luminosity is < 0. 1 , while only for the 
most extreme 20% it is between 0.5 and 0.7. Assuming the 
lower value, i.e. L s hcii/Jstar = 0.1, the total luminosities of 
the cocoon stars would be comparable to those of typical 
Galactic and Magellanic Clouds WCL stars, ~ 1O 55 J . 
But the lack of near IR spectral lines (especially in the 
J-band, see FMM99) suggests that Lshcii/Jstar must be 
large. If half the total stellar luminosity were re-cmittcd by 
the shell, and, for the sake of argument, if they were both 
Planckian with effective temperatures of 25,000 and 800 
K respectively, then the two energy distributions would 
intersect at ~ 1.37 fim, and the 1.25 /um continuum of the 
star would be ~ 3 x that of the shell. In this situation, the 
ratio ishcii/istar = f-0, already larger than for the most 
extreme DWCL stars in WHT. Increasing the stellar tem- 
perature to 35,000 K will move the intersection point to 
1.27 ^m, and the stellar continuum will be only 20% higher 
than the shell continuum. Under such circumstances the 
stellar emission line spectrum should still be easily visible 
in the J-band, unlike the results of FMM99. Increasing 
the shell luminosity to 10 x the (transmitted) stellar lu- 
minosity would be necessary for the shell continuum to 
dominate the stellar continuum at 1.25 /im. Thus, for the 
J-band spectrum to be dominated by the shell, the shell 
luminosity would have to be <; 90% of the total (shell + 
star), and the star has to be very hot. In this case, the 
UV optical thickness of the shells would have to be much 
greater than unity, which would probably be inconsistent 
with the optically thin shell in the near and mid IR, but 
it would bring the total luminosities roughly equal to the 
observed luminosities, and these values are rather low for 
typical DWCL stars. 

In summary, if the shells reprocess only a small part of 
the true stellar luminosity, then the true stellar luminosity 
is comparable to that of known DWCL stars, but the stel- 
lar spectrum should be observable in the J-band. If, on 
the other hand, the shells reprocess most of the luminos- 
ity of the central star, in which case our models are likely 
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to fail on optical thickness grounds, our models predict a 
total luminosity much lower than known DWCL stars. 

The physical size of the shells can be estimated from 
their luminosity and Tbb- Taking typical values of L = 
50, OOOLq and Tbb = 500K, we derive a physical diameter 
of 250 AU, which corresponds to and angular diameter of 
35 milliarcsec. This value should correspond roughly to the 
inner diameter of the shell. The outer diameter will be 50- 
100 x larger, or ~ 2-3". These values are very rough, but 
they are overall consistent with the observed extensions 
observed in the Palomar images, though the coincidence 
could be fortuitous. 

4-3. Other Interpretations 

Glass et al. (1999) suggested that the cocoon stars could 
be self-enshrouded massive stars of the type suggested by 
Bernasconi & Maeder (1996). Without discussing those 
models, we note that the measured shell luminosities in 
Table 4, if they represent the total stellar luminosity, cor- 
respond to stars of initial mass between 15 and 25M Q 
in the theoretical HR diagrams of Bernasconi & Maeder 
(1996). This would place them below the ~ 40M Q limit 
of the stars that would remain obscured on the Main Se- 
quence if those models were applicable. Secondly, if the 
cocoon stars are coeval with the rest of the cluster, their 
age of ~ 4Myrs would make them too old for this pro- 
posed evolutionary stage. 

Okuda et al. (1990) already discussed why other 
evolved stars are unlikely explanations of the nature of the 
cocoon stars. We add here that the luminosities we derived 
are somewhat too high for AGB stars, and also that the 
cluster is too young for AGB stars to have formed: stars 
much more massive than ~ 8M Q are still on the main 
sequence (FMM99). 

5. Conclusions 

New mid-infrared spectroscopic observations of the enig- 
matic cocoon stars in the Quintuplet Cluster were ob- 
tained in order to put further constraints on their nature. 
The high spectral/low spatial resolution spectroscopy did 
not reveal any feature that could be assigned to the cocoon 
stars themselves. The low spectral/high spatial resolution 
spectrophotometry was used to study the form of the sil- 
icate absorption in these stars, and we conclude that (i) 
the depth of the absorption is the same for all the stars, 
and (ii) there is no evidence for any silicate emission or 
absorption intrinsic to the cocoon stars. In fact, all of the 
absorption features are consistent with an origin in the 
molecular clouds and in the more tenuous ISM along the 
line of sight. 

We have examined carefully the hypothesis that the 
cocoon stars are extreme DWCL stars. The SEDs of the 
cocoon stars can be reproduced by the models of DWCL 
shells presented by WHT, but the high UV optical thick- 



ness needed to reprocess most of the luminosity of the 
central star suggests that these models might not be ap- 
plicable to such extreme stars. Furthermore, the luminosi- 
ties derived are significantly lower than those of known 
DWCL stars, and also these stars are extremely rare in 
the Galaxy, though they are concentrated in the GC. Yet, 
stars of the carbon sequence are necessary to produce fea- 
tureless spectra. 

The new observations do not provide any new con- 
straints in favour of any of the proposed identifications 
of the cocoon stars but they do provide new arguments 
against all three of the main hypotheses that have gener- 
ally been considered, namely dust-enshrouded youg stars, 
late-type (OH/IR or AGB) stars, or late-type WC stars. 
The nature of these stars remains an enigma. 
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